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Abstract 

Extensive studies on floral transition in nnodel species have revealed a network of regulatory interactions between proteins 
that transduce and integrate developmental and environmental signals to promote or inhibit the transition to flowering. 
Previous studies indicated FLOWERING PROMOTING FACTOR / {FPFl) gene was involved in the promotion of flowering, but 
the molecular mechanism was still unclear. Here, FPFl homologous sequences were screened from diploid Cossypium 
raimondii L. (D-genome, n = 1 3) and Gossypium arboreum L. genome (A-genome, n = 1 3) databases. Orthologous genes from 
the two species were compared, suggesting that distinctions at nucleic acid and amino acid levels were not equivalent 
because of codon degeneracy. Six FPFl homologous genes were identified from the cultivated allotetraploid Gossypium 
hirsutum L. (AD-genome, n = 26). Analysis of relative transcripts of the six genes in different tissues revealed that this gene 
family displayed strong tissue-specific expression. GhFPFI, encoding a 12.0-kDa protein (Accession No: KC832319) exerted 
more transcripts in floral apices of short-season cotton, hinting that it could be involved in floral regulation. Significantly 
activated APETALA 1 and suppressed FLOWERING LOCUS C expression were induced by over-expression of GhFPFI in the 
Arabidopsis Columbia-0 ecotype. In addition, transgenic Arabidopsis displayed a constitutive shade-avoiding phenotype that 
is characterized by long hypocotyls and petioles, reduced chlorophyll content, and early flowering. We propose that GhFPFI 
may be involved in flowering time control and shade-avoidance responses. 
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Introduction 

Cotton [Gossypium spp.) is one of the most important natural 
fiber crops in the world. In addition to its economic importance, 
cotton has attracted considerable scientific interest from plant 
breeders, taxonomists, developmental geneticists, and evolutionary 
biologists because of its unique reproductive developmental 
aspects and speciation history [1-3]. The genus Gossypium L. 
contains more than fifty species, which are cytogenetically 
differentiated into eight genomic groups (A-G, and K). Most of 
the species are diploid (n = 1 3), but five are allopolyploids (n = 26), 
originating from an interspecific hybridization event between A- 
and D-genome diploid species. Not only two diploids of A- 
genome, Gossypium herbaceum L. and G. arboreum L., but also two 
allopolyploids of AD-genome, G. hirsutum L. and G. barbadense L. 
were domesticated by humans for their fiber demands [3] . Because 
of environmental pressures, such as land use and climatic change, 
the earliness of cotton has become a vital subject for plant 
breeders. Several traits co-operate to influence the early ripeness of 
upland cotton (G. hirsutum L.), with flowering time being especially 
important. In the seed crop, floral transition is a key developmen- 
tal switch in the life cycle of cotton as it contributes to the 



production of dry matter. Furthermore, shifting of the seasonal 
timing of reproduction is a major goal of plant breeding research 
as it wiU produce novel varieties better adapted to local 
environments and effects of climate change [4] . 

Plant growth originates from a small number of undifferen- 
tiated cells called meristems. The apical meristems being 
indeterminate or determinate contribute to the fate of shoot 
architecture. Indeterminate apical meristems retain a population 
of vegetative stem cells indefinitely which perform tissue and 
organ differentiation below and on the flanks of the main-stem. 
Meanwhile the determinate apical meristems undergo terminal 
differentiation, commonly in a flower or inflorescence. In annual 
plants, the floral induction process occurs when vegetative shoot 
meristems develop into inflorescence meristems, and give rise to 
flowers [5] . The use of the model species Arabidopsis thaliana and 
Antirrhinum majus has led to significant progress in the 
understanding of the floral transition. Numerous genes involved 
in the control of flowering time have been identified, and the 
roles they played in molecular and genetic pathways were also 
characterized. Previous studies uncovered that four main floral 
pathways, namely vernalization, photoperiod, gibbereUin, and 
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autonomous pathways converged to regulate floral integrator 
genes such as LEAPT (LFT), APETALA 1 [API), FLOWERING 
LOCUS T {FT) and SUPPRESSOR OF OVEREXPRESSION OF 
CONSTANSl {SOCT), which in turn could activate genes 
required for reproductive development [6]. Another transcrip- 
tion factor FLOWERING LOCUS C {FLQ containing MADS 
domain was an inhibitor of flowering, acting as an important 
convergence point for the autonomous and vernalization 
pathways in Arabidopsis thaliana [7-9]. 

FLOWERING PROMOTING FACTOR 1 {FPFl) gene was 
originally understood on account of its role in flowering. Over- 
expression of FPFl (Y 11988) led to early flowering in Arabidopsis 
[10]. FPFl is expressed in apical meristems immediately after the 
photoperiodic induction of flowering in long-day plants that could 
flower in response to long days [10]. Previous studies indicated 
that FPFl might pla)' an important role in modulating the 
competence of apical m(;rist(-ms to rapidly respond to the floral 
meristem identity genes APland LFT. During the transition to 
flowering in Arabidopsis, FPFl is normally activated at a similar 
time as LFT, and earlier than API. Over-expression of AtAPl and 
AtFPFl shows a synergistic effect in the shortening of the time to 
flowering both under long-day and short-day conditions [1 1]. Two 
closely related genes, FLPl and FLP2 (FPFl-Uke genes), have been 
identified in Arabidopsis. Constitutive over-expression of each gene 
causes earlier flowering under both long and short day conditions 
[12]. Up till now, homologous genes of FPFl have also been 
characterized in rice {Otyza sativa) and tobacco {Nicotiana tabacum) 
[13,14]. OsRAAl, a homolog of FPFl in rice, shares similarity of 
58% with AtFPFl at amino acids level. Evidence revealed that 
over-expression of OsRAAl causes pleiotropic pheno types in 
transgenic rice plants, including altered leaf shape, heading time 
and root development [13]. 

In summarize, AtFPFl / OsRAAl gene family takes parts in 
several aspects of plant development. Currently, litde is known 
about the underlying mechanism of these, or of any crossover 
between this gene family and other floral molecules. This paper 
described our work to dissect these mechanisms. In our previous 
study, a high-quality, normalized, full-length cDNA library with 
a total of 14,373 unique ESTs was generated to provide 
sequence information for gene discovery related to flower 
development in upland cotton [15]. The publication of G. 
raimondii L. genome sequence by the Cotton Genome Project in 
2012 facilitates gene excavation vital for the genetic improve- 
ment of cotton quality and productivity, as well as serving as a 
reference for the assembly of the tetraploid G. hirsutum L. 
genome [16]. G. raimondii L. genome (http://cgp.genomics.org. 
cn/) and G. arboreum L. genome databases (unpublished) were 
screened to identify and compare homologs of AtFPFl from the 
two diploid cotton species. The genes were cloned and 
characterized from G. hirsutum L. as one of the major cultivated 
species. In addition, their expression pattern in cotton, and 
ectopic expression in Arabidopsis were analyzed. 

Materials and Methods 

Plant Material and Growth Conditions 

A genetic standard line TM- 1 (G. hirsutum L.) and a short-season 
cotton \'ariety CGRI 36 (G. hirsutum L.) were grown in a 
greenhouse (Anyang, China) with an optimal temperature of 
28°C. In order to minimize experimental error, two varieties were 
planted in the soil possessing equal fertility and managed using 
standard agricultural practices. Shoot apices from about three 
hundred plants of the two varieties were harvested when three true 
leaves came out. Roots, stems, young leaves, mature leaves. 



flowers, and fibers from mature CCRI 36 and TM-1 plants were 
collected and mixed together. 

Seeds of Arabidopsis thaliana (Columbia-0 ecotype: wild-type or 
GhFPFl over-expressing transgenic plants) were sterilized in 
75% (v/v) absolute alcohol for 30 s and \0% Yf^O^ (v/v) for 
10 min. After four washes in sterihzed double-distilled water 
(ddHgO), seeds were sown onto sterilized 1/2 MS sohd medium 
(PH 5.8) containing 1.5% sucrose, 0.7% agar (and 50 mg/L 
kanamycin for transgenic lines). 10-cm petri plates containing 
medium and seeds were chilled for 48—72 h at 4°C in the dark, 
and then transferred into an illuminating incubator with 
100 |a,mol-m~^ s~' fluorescent light at 22°C. Two weeks later, 
normal seedlings were transplanted into the soil in a growth 
chamber which provided the plants with continuous 
150 |J,mol-m ^ s ' fluorescent light at room temperature 
22°C. As for the length of h^'pocotyl and petiole, chlorophyll 
content and flowering time assay, seeds were directly sowed into 
the soil in the growth chamber. Ten days later, seedlings were 
transplanted into flowerpots grown under the same conditions. 
AH Arabidopsis plants were grown under long-day conditions with 
high red to far red (R/FR ratio: 4.5) provided by fluorescent 
lamps. 

Gene Screening and Sequence Comparative Analysis 

Candidate sequences were obtained by performing Blast 
searches of genome databases of G. raimondii L. (http://cgp. 
genomics.org.cn/) and G. arboreum L. (unpubhshed) using Arabi- 
dopsis thaliana FPFl gene as a query sequen[:e. Multiple sequence 
alignment of the deduced amino acids was performed using the 
software ClustalX2 (http://www.ebi.ac.uk). The phylogenetic tree 
was constructed by Molecular Evolutionary Genetics Analysis 
(MEGA) software 4.1 (http://www.megasoftware.net). Sequence 
of promoter was analyzed in the database of plant cis-acting 
regulatory DNA elements (PlantCARE) (http://bioinformatics. 
psb.ugent.be/ webtools/ plantcare/html/). 

Gene Isolation and RNA Analysis 

Different organs from upland cotton and Arabidopsis were 
harvested, and immediately plunged into liquid nitrogen. Total 
RNA was isolated from cotton and Arabidopsis using the EASY spin 
plus RNA reagent kit RN38 (AIDLAB, Beijing, China) following 
manufacturer's instructions. Poly (dT) cDNA was prepared from 
the total RNA using the Superscript III First-Strand Synthesis 
System (Invitrogen, USA) according to the manufacturer's 
protocol. Genomic DNA was extracted as described [17]. 

Open reading frames (ORFs) of GhFPFland homologous genes 
were amplified from G. hirsutum L. CCRI 36 using the primers 
described in Table SI. The running conditions were as follows: 
denaturation at 94°C for 2 min, followed by 33 cycles at 94°C for 
30 s, annealing at 58°C for 1 min and extension at 72°C for 30 s. 
To obtain the 5 '-terminal and 3 '-terminal sequences of GhFPFl, 
5'- and 3'-RACE were performed using the SMART RACE 
cDNA Amplification Kit (Clontech, USA) according to the Kit's 
protocol in conjunction with primers GPSl and GPS2 (Table SI). 

Quantitative polymerase chain reactions (QRT -PCRs) were 
performed to measure relative expression of homologous genes 
using specific primers (Table S2), with Histone 3 (AF024716) to be 
an internal control [18]. For Arabidopsis, AtUBQ5 (AT3G62250) 
was used as the control. Transcriptional changes of flowering 
related genes AtFPFl, AtLFT, AAPl, AtCO, AtFT, AtSOCl, AtFLC 
and AtPHTB were examined in the wild as well as transgenic 
plants. Each qRT-PCR experimental condition was independentiy 
repeated three times and in each of these three biological 
repetitions, three technical replicates were made. AH of the 
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Table 1. Pair-wise alignment between G. raimondii L. and G. arboreum L. orthologous sequences on amino acid and nucleic acid 
levels. 



G. raimondii L DF10021325 DF10039615 DF1000745S DF10000980 DF10029023 DF10009151 

G. arboreum L Garb08271 Garbl7683 Garbl9766 Garb07734 Garb34130 Garb22630 

Similarities of nucleic acid/amino acids (Identity %) 97/100 99/100 95/98 99/99 97/99 84/86 



doi:l 0.1 371/journal.pone.0091 869.t001 

amplification of interested genes was analyzed on ABI 7500 system 
(Applied Biosystems, USA) with SYBR Green I (with Rox) 
reagents to detect the target sequences. The running conditions 
were as follows: holding stage at 50°C for 2 min, 94°C for 10 min, 
followed by 40 cycles at 95°C for 15 s, 60°C for 1 min. Then a 
melting curve was performed from 65 to 95°C to verify the 
specificity of the amplified product. QRT-PCR data was processed 
to measure the relative expression of genes in accordance with the 
2 method [19]. 

Jasmonic Acid and Salicylic Acid Treatment 

Jasmonic acid (JA; Sigma, USA) was dissolved in a small 
quantity of ethanol, and further diluted with ddHaO to a 
concentration of 200 |J.mol-L Salicylic acid (SA; Sigma, USA) 
was directly dissolved into heated ddHgO water to a final 
concentration of 400 |.tmol-L '. Upland cotton cultivar CCRI 
36 was planted in the greenhouse with an optimal temperature of 
28°C. One-month-old plants were sprayed with JA and SA 
solutions when four true leaves had developed [20,21,36]. The top 
two leaves were removed at 0, 0.5, 6, 12, 24, and 48 h time-point. 
Harvested leaves were immediately plunged into liquid nitrogen, 
and stored at — 80°C for later RNA extraction and relative 
transcript analysis. Nonexpressor of Pathogenesis-Related Genes 1 
{GhJVPRl, DQ409173), a SA- and JA-inducible gene in upland 
cotton was used as a positive control gene here [21]. 

Vector Construction and Genetic Transformation 

The open reading frames of GhFPFl (GenBank Accession No. 
KC8323 1 9) was cloned into the binary vector pBI 1 2 1 using an In- 
Fusions ™ Advantage PGR Gloning Kit (Glontech, USA) via 
BamHl and Sad sites (New England BioLabs, USA). The 
recombinant plasmid containing GaMV35S::GAi^P-F7 was intro- 
duced into Agrobacterium tumefaciens (strain LBA4404) and then 
transformed into Arabidopsis (Columbia-0 ecotype) according to the 
method of floral dip [22]. Transformants were selected on 1/2 MS 
medium containing 50 mg 1 ' kanamycin, and further confirmed 
at both the genomic DNA and transcriptional mRNA level. To 
confirm that GhFPFlhad been integrated into the Arabidopsis 
genome, RT-PCR was performed with the 35S promoter primer 
(forward) and the GhFPFl -speciSic primer (reverse). To detect the 
mRNA expression oi GhFPFlin T3 transgenic lines COL-3, GOL- 
4 and GOL-7, 18-day-old intact plants were sampled for 
quantitative RT-PGR. 

Flowering Time Calculation and Chlorophyll Content 
Determination 

At least twenty individual plants of the homozygous T3 
transgenic lines and wild-type, were governed under the long- 
day conditions (16 h light/8 h dark). Flowering time was 
measured by counting the numbers of rosette and cauline leaves 
and days of the appearance of the first flower. Statistical 
significance analysis was conducted using the software Sigma Stat 
3.5 (Systat Software, San Jose, GA, USA) by One-Way ANOVA 



analysis. Leaves to determine the content of chlorophyll were 
harvested from 24-day-old Arabidopsis. Ghlorophyll was extracted 
through 95% acetone-alcohol according to the method as 
described previously [23] and monitored in a 96-weU microplate 
reader (BioTek, USA). 

Ethics Statement 

We did not make use of human or vertebrate animal subjects 
and/or tissue in our research. 

Results 

Gene Screening and Sequence Comparative Analysis 

To investigate the homologs of FLOWERING PROMOTING 
FACTOR 1 [FPFl] in cotton, the genome databases of G. raimondii 
L. and G. arboreum L. were screened with the coding region of 
Arabidopsis FPFl gene as the reference sequence. She sequences, 
which covered complete ORFs, were selected from each database. 
The twelve ORFs were predicted to encode small proteins of 99 to 
1 1 3 amino acids in length. Pair-wise alignment of these predicted 
proteins suggested that they shared 50-68% similarity to the 
AtFPFl protein. In addition, these proteins were rich in the three 
amino acids Ser, Val, and Leu. 

Further, pair-wise alignment of the amino and nucleic acids was 
conducted to explore the differences between G. raimondii L. and G. 
arboreum L. orthologous genes. Nucleic acid alignment revealed 
that five of the six orthologous pairs exhibited 95-99% similarity 
to each other, and the other pair DF10009151- Garb22630, had 
84% identity (Table 1). Two pairs among them, DF10021325- 
Garb08271 and DF10039615-Garbl7683, owned the same 
deduced amino acids sequence because of codon degeneracy. 
One or two amino acids were different between the three pairs of 
orthologous DF10000980-Garb07734, DF10029023-Garb34130 
and DF10007455-Garbl9766. The remaining pair, DF10009151- 
Garb22630, showed greater divergences, with just 86% homology 
at amino acid level. Comparative analysis revealed that a short 
length of nucleic acids was missing from Garb22630 (Figure SI). 

Isolation of FPFl Homologous Genes in G. hirsutum L. 
and Sequence Analysis 

Six genes were identified from GGRI36 and named as GhFPFl, 
GhFLP-1, GhFLP-2, GhFLP-3, GhFLP-4, and GhFLP-5 correspond- 
ing to the sequences of DF10009151, DF10007455, DF10029023, 
DF10021325, DF10000980, and DF10039615 in G. raimondii L., 
respectively. The small FPFl gene family has been characterized 
in several species including Arabidopsis thaliana, white mustard 
[Sinapis alba), rice [Oryza saliva), tobacco (Mcotiana tabacum) and 
maize {^ea mays). Multiple alignment (Figure lA) of amino acid 
sequences revealed that distinctions between the homologous 
proteins occurred in the N-terminus relative to C-terminus. 
Results suggested that there were three conserved domains present 
in the protein family. The first motif, -LGWERY- was located in 
the middle section of the protein, while the second and third 
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conserved motifs, -D/HLISLP- and -MY/FDIWKN-, were 
found closer to the C-terminus. Phylogenetic analysis of the 
FPFl gene family indicated that these proteins could be divided 
into three clades, represented by NtFPFl, ^niFPFl / OsRAAl and 
AtFPFl. GhFPFl, GhFLP-1, GhFLP~2, and GhFLP-3 were placed 
into JVtFPFl clade; GhFLP-4 and GhFLP-5 were in the same clade 
with AtFPFl. None of them were placed into the Z"iFPFl / OsRAAl 
branch (Figure IB). 



GhFPFl had Higher Transcriptional Levels in the Floral 
Apices of CCRI 36 

QRT-PCR was used to profile relative expression of above six 
genes in different tissues of upland cotton. Roots, stems, leaves, 
flowers and fibers are mixed samples from CCRI 36 and TM-1 
while floral apices T and C represent floral apices from CCRI 36 
and TM-1, respectively (Figure 2). GhFPFl gene family displayed 
tissue-specific expression because abundant transcripts of the six 
genes were found in roots, floral apices, flowers, and stems, but 
were barely detectable in leaves or fibers. GhFLP-3, GhFLP'4 and 
GhFLP-5, in particular, had only some but not much expression in 
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Figure 1. Multiple sequence alignment of FPF1 protein family from G. hirsutum L., and other species. A. Multiple alignment of FPF1 
protein sequences in several species. AtFPFl (Yl 1988) and AtFLPl (AL353995) are Arabidopsis thaliana genes; SaFPFl (Y1 1987), NtFPFl (AY496934), 
ZmFPFI (ACG44143) and OsRAAl (AY659938) are from Sinapis alba, Nicotiana tabacum, Zea mays and Oryza sativa. B. Phylogenetic tree of the FPFl 
proteins in the above plants as determined by the IVIEGA 4.1 software package. 
doi:1 0.1 371/journal.pone.0091 869.g001 
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floral apices relative to other tissues. More importantly, we focused 
on the contrastive analysis of gene expression in floral apices of 
CCRI 36 (a short-season cotton variety) and TM-1 (a genetic 
standard line). Results uncovered that GhFPFl had more than 
four-fold transcript levels in CCRI 36 than in TM-1. Higher 
expression of GhFPFl in the short-season cotton ,sugg(-sted that it 
was the most possible FPFl orthologous gene as AtFPFl involved 
in the promotion of flowering. 

Gene Structure Analysis of GhFPFl and its Response to JA 
and SA 

A b'- and 3'- RACE strategy was performed to gain 
transcription initiation and termination sites of GhFPFl. A full- 
length cDNA of 701 bp composed of 56 bp 5'-UTR, 315 bp 3'- 
UTR and 330 bp ORF was isolated from the cDNA pool of one- 
week-old seedlings (Figure 3A). Comparison of genomic and 
cDNA sequence revealed that there was no intron. Sequences of 
GhFPFl and GhFPFl -like genes were submitted to NCBI (Acces- 
sion number: KC832319, KF830866-KF830870). 

A 2076 bp promoter region of GhFPFl was identified from the 
genomic DNA of upland cotton, containing one and three 
response elements of methyl jasmonate and salicylic acid 
respectively. Treatment assay and the following qRT-PCR were 
performed to evaluate whether GhFPFl could be regulated by JA 
and SA plant hormones (Figure 3B). NPRl [nonexpressor of 
pathogenesis-related genes 1), a pathogen-related gene regulating SA- 
dependent defense response, systemic acquired resistance, and 
mediating crosstalk between SA and JA [24,25], was used as a SA- 
and JA- inducible reference gene. According to figure 3B, GhFPFl 
and GhNPRl exhibited the similar expression trend when plants 
were treated with SA andJA, respectively. After 12h treating with 
SA, the expression of GhNPRl reached its peak but peak 



expression of GhFPFl occurred 12h later. However, GhFPFl was 
responsive to JA earlier than GhNPRl. GhFPFl transcripts were 
decreased 0.5h later and began to rise till 12h after the JA 

treatment, but the transition points of expression trend of GhNPRl 
was lagging. GhFPFl showed different expression fluctuations 
exposed to exogenous SA and JA but it was positively regulated by 
the two phytohormones over most time of 48h, suggesting that 
GhFPFl could be responsive to both JA and SA, the latter later. 

Over-expression of GhFPFl Promoted Flowering in 
Arabidopsis 

To discuss whether GhFPFl could modulate flowering time in 
plants, the open reading frame of GhFPFlwAti transformed into the 
Arabidopsis Columbia-0 ecotype under the control of the cauliflow- 
er mosaic virus 35S promoter. Seven independent transgenic lines 
were obtained, and the GhFPFl transcript was analyzed in three 
transgenic lines COL-3, COL-4, and COL-7 (Figure 4D). 
Homozygous T3 lines and the wild-type were grown under 
long-day conditions of 16 h light/8 h dark. Twenty-four days after 
sowing, most plants of the transgenic lines had started bolting, 
while wild-type plants remained in vegetative growth with six 
rosette leaves developed (Figure 4A). Six days later, the transgenic 
plants flowered in succession, whereas the wild-t)'pe plants had just 
completed their basal r()S('tte dcvc'lopmciit at this point (Figure 4B 
and C). To determine flowering time, the number of rosette and 
cauline leaves and days of the first flower opening up were counted 
(Table 2). Results turned out that transgenic plant produced flower 
buds six days earlier than 36.8 days to flowering in the wild-type 
on average. This corn^sponds to a reduction in the number of 
leaves from 15.4 in wild-type to the fewest 1 1.5 in the transgenic 
line. Data demonstrated that over-expression of GhFPFl in 
Arabidopsis promoted flowering under inductive photoperiods. 
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Figure 2. Expression patterns of GhFPFl, GhFLP-1, GhFLP-2, GhFLP-3, GhFLP-4, and GhFLP-Sm G. hirsutum L. Relative expression of GhFPFl 
and its homologs were measured in different tissues of upland cotton CCRI 36 (a short-season cotton variety) and TIVI-1 (a genetic standard line) using 
qRT-PCR. Roots, stems, leaves, flowers and fibers stand for mixed samples from CCRI 36 and TIVI-1. Floral apices from CCRI 36 and TM-1 were 
harvested and named as floral apices T and C respectively. Error bars represent standard deviation (SD). 
doi:1 0.1 371/journal.pone.0091 869.g002 
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Figure 3. Gene structure of GhFPFI, and its response to plant hormones. A. The gene structure of GhFPFI. JA and SA represent response 
elements of JA and SA; TSS, transcriptional start site. B. GhFPFI and GhNPRI expression profiles in the first forty-eight hours after treatment with SA 
and JA. GhNPRI (DQ409173), a known JA- and SA-inducible gene, was used as a positive control here. Error bars represent SD. 
doi:1 0.1 371 /journal.pone.0091 869.g003 



To gain further insight into how GhFPFI regulates the floral 
transition, we examined the expression levels of six genes 
associated with the promotion of flowering of Arabidopsis: FPFl, 
LFT, API, CONSTAMS {CO), FT, SOCl, and the flowering 
repressor, FLC both in wild-type and transgenic plants of 
vegetative (18-day-old) and bolting on (24-day-old). Relative 



qRT-PCR indicated that the level of AtFLC transcripts in 
transgenic plants of 18-day-old was one-eighth of that in wild- 
type plants under long-day conditions. Meanwhile, AtSOCl and 
AtAPl transcripts were up-regulated by two to three-fold relative to 
the wild-type. There was no obvious change in the relative 
transcript amount oiAtFPFl, AtLFF, AtFT, or AtCO (Figure 5A). In 




35S::GhFPF1 



Figure 4. over-expression of GhFPFI promoted flowering in Arabidopsis. A. 24-day-old plants grown under long-day conditions. The wild- 
type control developed six small rosette leaves (left) by the time that the GhFPFI transgenic plant had started bolting on (right). B and C. One-month- 
old plants grown under long-day conditions. The wild-type control was still at the vegetative stage (left), whereas the transgenic plant(s) had 
flowered (right). D. Relative transcriptional analysis of GhFPFI in 18-day-old transgenic Arabidopsis lines COL-3, COL-4 and COL-7. 
doi:1 0.1 371/journal.pone.0091 869.g004 
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Table 2. Flowering time under long-day conditions, as measured by days to flowering and number of basal rosette and cauline 
leaves. 



Line NO 


Flowering time (days after sowing) 


P value 


Rosette leaves + cauline leaves 


P value 


WT 


36.8±0.47 




15.4±0.55 




COL1 


30.0±0.44* 


<0.05 


11.5±0.57* 


<0.05 


COL2 


30.9±0.34* 


<0.05 


12.4±0.42* 


<0.05 


COL3 


30.3±0.53* 


<0.05 


11.9±0.42* 


<0.05 


COL4 


33.3±074* 


<0.05 


12.7±0.41* 


<0.05 


COLS 


32.2±072* 


<0.05 


13.3±0.57* 


<0.05 


COL6 


32.5±0.67» 


<0.05 


13.3±0.56* 


<0.05 


COL7 


30.7±0.66* 


<0.05 


12.2±0.32» 


<0.05 



Asterisks indicate significant variation differences between the wild-type and eacfi 35S::GhFPFl transgenic population line as determined by one-way ANOVA analysis. 
Means with SD from twenty plants of every line were shown. 
doi:1 0.1 371 /journal.pone.0091 869.t002 



genetic divergence might occur during speciation for five genes 
liad higlier synonymous changes between the two species. 

Previous studies have indicated that for each gene studied, 
allotetraploid species of cotton should liave two homologs which 
represent descendants from the A-genome, and D-genome donors, 
evolving independently at the time of polyploidy formation [26] . 
Meanwhile, on allopolyploidy, genomic changes will take place, 
including chromosomal rearrangement and changes in gene 
expression. During the cloning and identification of the FPFl 
homologous genes in G. hirsutum L., at least eight single clones were 
sequenced for every gene. It was discovered that five genes 
contained cDNA sequences identical to D-genome sequences from 
G. raimondii L., whereas only one gene GhFLP3 was found to have 
the same sequence as the A-genome sequence from G. arboreum L, 
The finding revealed that this gene family exhibited subgenome- 
specific expression bias to D-subgenome. 

When plants undergo the transition to flowering, the vegetative 
shoot apical meristem will be transformed into an inflorescence 
meristem. Inflorescence meristems can respond to both environ- 
mental and endogenous flowering signals to give rise to floral 
meristems, which go on to produce the various floral organs in 
succession [27]. Axillary buds of upland cotton can be either 
vegetative or floral buds which develop into vegetative shoots or 
reproductive shoots, respectively. Two types of axillary buds are 
distinct from each other for apical meristems of the vegetative buds 
are small and either domed or tapered in shape, with one or two 
layers of tunica cells but the apical meristems of the floral buds are 
large and columnar, with two or three layers of tunica cells in floral 
buds with flat surfaces [28] . It was reported that the short-season 
cotton varieties of CCRI16 and CCRI36 initiated morphological 
dififerentiation of the floral bud at the developmental stage of two 
true leaves flattened, whereas the late-maturing CCRI12 activated 
this process when there were three true leaves expanded [28,29]. 
TM-1, a genetic standard line as well as a late-maturing variety 
was also found to begin morphological differentiation of the floral 
bud when three true leaves were expanded (Figure S2). To select 
specific genes that are involved in the floral development of short- 
season cotton, the relative expression of GhFPFl and GkFPFl -like 
genes in the floral apices of CCRI 36 and TM-1 were analyzed 
and compared at the developmental stage of three true leaves 
expanded. Transcripts of GhFPFl were discovered to be more 
abundant in the floral apices of CCRI 36, suggesting that GhFPFl 
might be involved in the floral regulation of short-season cotton. 



24-day-old transgenic plants collected at the same time of 1 1 am (3 
hours exposure to light), the transcripts of AtFPFl, AtLFT, AtFT, 
AtCO were up-regulated by four times more or less, but AtSOCl 
transcripts remained at the same levels as in 18 days (Figure 5A 
and B). Remarkably, AtAPl was activated highly with an increase 
of 18.5-fold compared to the wild-type; meanwhile, the expression 
of the flowering repressor AtFLCwas almost completely suppressed 
due to the ectopic expression of GhFPFl. Combining these, it 
could be hypothesized that early flowering conferred by GhFPFl 
over-expression in Amhidopsis might be mediated through AtAPl 
and AtFLC. 

Over-expression of GhFPFl Triggered Shade Avoidance 
Syndrome (SAS) in Arabidopsis 

Except accelerated flowering time, compared with wild-type, 
transgenic Arabidopsis generated longer hypocotyls and petioles 
(Table 3), as well as the upward of movement of leaves. Also 
chlorophyll content was found to be reduced in transgenic plants 
(Figure 6). Taken together, transgenic plants were recognized as 
so-called shadow avoidance syndrome (SAS). The elongated 
appearance and early flowering response were similar to 
phenotype of phjiB mutants. PHYB, acting as the major 
phytochrome in light-grown plants played a predominant role in 
shade avoidance syndrome [48]. Since the plants were provided 
sufficient fluorescent light with red to far-red ratio 4.5, transcripts 
of PHYB were measured using qPCR in transgenic and wild 
Arabidopsis further. The result revealed that expression of PHYB in 
transgenic plants was decreased by more than fifty percent 
(Figure 6E). 

Discussion 

GhFPFl gene belongs to a novel gene family that seems to be 
conserved in both higher and lower plants. Until now, as were 
characterized in several species, members of this gene family were 
short in length as well as lacking in intron of their genomic 
sequences [10-12]. Twelve FPFl homologs were identified from 
the diploid cotton genomic databases of G. raimondii L. and G. 
arboretum L. Orthologous sequences from the two cotton species 
were compared with each other, suggesting that nucleic acid 
sequences of the six pairs of orthologs were distinct, though two 
pairs of orthologous genes possessed the same deduced protein 
sequence as a result of codon degeneracy. High constraints of 
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-6 ■> 



AtFPFI AtLFY AtCO AtFT AtAPI AtSOCI AtFLC 



Figure 5. Relative qRT-PCR analysis of genes involved in floral 
regulation in Arabidopsis. The AtUBQS gene was used for calibration 
and the histogram was drawn based on the log2 scale of the ratio of 
gene expression in transgenic plants relative to wild-type. Entire 18-day- 
old plants of vegetative (A) and the above-ground parts of 24-day-old 
plants of starting bolting on (B) were analyzed. Interested genes were 
AtFPFI (AT5G24860), AtLFY (AT5G61850), AtCO (AT5G15840), AtFT 
(AT1G65480), AtAPI (AT1G69120), AtSOCt (AT2G45660), and AtFLC 
(AT5G10140). All plants were grown under long-day conditions of 16 h 
light/8 h dark. Data in graph were mean values with standard deviation 
(error bar) from three replicates. 
doi:1 0.1 371 /journal.pone.0091 869.g005 

In the study we found that GhFPFl was predominantly 
expressed in roots and floral apices, the former more. In Arabidopsis 
thaliana, a number of genetic pathways controlling flowering time 



have been identified. In fact, most of the genes involved were 
preferentially expressed in the shoot apical meristem (SAM) and 
the root tip, such as PHTA, CRT2, FLC and so on, but, surprisingly, 
only a few were expressed preferentially or exclusively in leaves for 
example FT [7,30,31]. Day length and light quality are essentially 
perceived by photosynthetic organs, expanded leaves and stem, 
whereas water and mineral availability are perceived by the roots. 
The root system is presumably capable of reacting to the critical 
environmental changes and, as a result, influences the flowering 
process to some extent. In Sinapis, analyses of changes in the 
contents of phloem and xylem saps during the floral transition 
have disclosed a complex shoot-to-root-to-shoot signaUing loop 
involving both nutrients and hormones [32,33]. Higher transcripts 
of GhFPFl in roots implicates that it could be involved in multiple 
functions in Gossypium hirsutum L. Moreover, members of GhFPFl 
gene family displayed tissue specific expression, which is not rare. 
Most genes controlling flowering time were expressed across a 
wide range of organs and tissues, but a survey of available data on 
their spatial expression patterns revealed that many genes showed 
preferential expression in more limited areas such as SAM and 
RAM [34]. 

Jasmonic and salicylic acid response elements were found in die 
promoter region of GhFPFl, which suggested that GhFPFl may be 
regulated by the plant hormones JA and SA. Response of GhFPFl 
to JA (wound signal) and S A treatments could hint involvement of 
GhFPFl in plant defense responses. Previous studies implicated 
both positive and negative roles of SA in affecting flowering time 
via influencing the expression of flowering regulatory genes FLC 
and FT [35-38] . Also SA can links flowering time with some stress 
signalling coming from defense responses or poor-nutrition 
[36,38]. Recent studies provided new insights into the mechanisms 
of salicylic acid (SA) perception and NPRl was proposed to be with 
.NPR3/.NPR4, resembling the multi-receptor of SA in diverse 
immune responses such as basal defense, systemic acquired 
resistance establishment, and effector-triggered immunity (ETI) 
[39]. GhFPFl and GhJVPRl shared the similar expression trend to 
exogenous SA but the association between FPFl, SA, plant 
defense responses and flowering time requires further investiga- 
tion. Little has been reported about the relationship between JA 
and the regulation of flowering time. Thus, it seems that JA may 
not have a direct effect upon the transition to flowering. Because of 
domestication, upland cotton has become a compact day-neutral 
and an annual row-crop from a lanky photoperiodic and perennial 
plant. Moreover some growth characteristics of cotton are 
diflFerent from other species as vegetative growth will proceed 
after the initiation of reproductive growth. Also flowering and fruit 
set are not synchronized but continue through the growing season. 
These competing sinks in upland cotton may give rise to flowering 
mechanism different from Arabidopsis. 

Ectopic expression of GhFPFl in Arabidopsis led to earlier 
flowering, and a decreased number of rosette and cauline leaves. 



Table 3. Statistics of hypocotyl and petiole lengths of wild and transgenic plants grown under long-day conditions with 
fluorescent lamps. 





Line NO 


Hypocotyl length (mm) 


P value 


Petiole length (mm) 


P value 


WT 


3.8±0.04 




7.0±0.18 




35S::GhFPF7 


6.5 ±0.07" 


<0.01 


11.2±0.21** 


<0.01 



Hypocotyl lengths of wild and transgenic lines were quantified from thirty ten-day-old plants. Petiole lengths of true leaves of the basal rosette were measured and the 
data were gathered from thirty wild type and transgenic 24-day-old plants, respectively. Mean values (±SD) were shown and statistical analysis was evaluated by the 
same way as table 2. 
doi:l 0.1 371/journal.pone.0091 869.t003 



PLOS ONE I www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e91869 



Characteristics of FPFl Ortliologs in Gossypium L. 





D 




WT GhFPFI 



Figure 6. Over-expression of GhFPFI led to sKiade-avoidance responses in transgenic plants. A. Phenotype of hypocotyl of ten-day-old 
wild type (left) and transgenic plants (right) grown under long-day conditions. B. The seventh rosette leaves of wild type (left) and transgenic (right) 
24-day-old plants grown under the same conditions as (A). C. Shade-avoidance responses in 24-day-old transgenic plants (right) which grew fast with 
upward leaves. All Arabidopsis plants were grown under long-day conditions with high red/far red (R/FR ratio: 4.5) light provided by fluorescent 
lamps. D. Chlorophyll content of leaves in transgenic plants was lower than that in wild-type and the difference was very significant (P<0.01) 
assessed by T-test. E. The transcript levels of AtPHYB (AT2G18790) in wild-type and GhFPFI over-expression transgenic plants. The AtUBQS gene was 
used as calibrator. 

doi:1 0.1 371/journal.pone.0091 869.g006 



When compared to the wild-type, transgenic plants had an 
increased expression of AtAPl, and suppressed expression of 
AtFLC. The MADS-domain transcription factor API, acts as a 
floral meristem identity gene that controls the onset of Arabidopsis 
flower development. API expression is first observed throughout 
the emerging floral primordia, and is later confmed to the outer 
whorls of floral buds, where it is also involved in the specification 
of sepals and petals [40,41]. Another MADS-box transcription 
factor, FLC is a major repressor of flowering in Arabidopsis. It binds 
to the first intron oi FT, and the promoter of SOCl, in each case 
inhibiting transcriptional activity. The FT protein interacts with 
FD to stimulate the activity of API. SOCl can bind to the 
promoter of LFT to activate its transcription. The actions of API 
and LFY promote the development of the inflorescence meristem, 
which leads to the production of flowers [42-44]. Ectopic 
expression of GhFPFI in 18-day-old Arabidopsis caused slight 
increases in AtAPl and yli50C7 expression levels. AtAPl transcript 
was induced to a much higher level six days later, when the 
transgenic Arabidopsis had started bolting. ^tfOCiexpression 
remained at the same levels as before. AtFLC expression was 
suppressed at an extremely low level at both time points. This 
reveals that the promotion flowering of over-expression of GhFPFI 
in Arabidopsis is possibly dependent on AtAPl and AtFLC. Previous 
studies had confirmed that the FPFl gene family took a positive 
effect on flowering time regulation [10,13,14]. Until now, little was 
known about the molecular mechanism of FPFl in floral 
regulation pathways. It is expected that the role of FPFl in the 
transition to flowering will be the focus of future research. 



To grow and develop optimally, all organisms need to perceive 
and process information from their environment. As sessile 
organisms, plants need to sense and respond to external stimuli 
more than most organisms. Therefore, plants have to adapt their 
developmental pattern to the environmental changes to ensure 
survival and reproduction. Light influences every developmental 
transition from seed germination and seedling emergence to 
flowering. For shade-intolerant plants, such as Arabidopsis thaliana, a 
reduction in the red to far-red (R: FR) ratio of incoming radiation, 
which is caused by absorption of red light and reflection of far-red 
radiation by canopy leaves, signals the proximity of neighboring 
plants and triggers the shade avoidance syndrome (SAS) [45]. A 
common phenotype of the SAS is re-allocation of energy resources 
from storage organs to stems and petioles so that the plant 
outgrows its competitors. Other responses induced by reduction in 
R: FR ratio include increased leaf angle, accelerated leaf 
senescence and reduced deposition of fixed carbon to storage 
organs [46] . Shade avoidance in higher plants is regulated by the 
action of multiple phytochrome (phy) species that detect changes 
in the red to far-red ratio (R: FR) of incident light to initiate a 
redirection of growth and an acceleration of flowering [47]. 
Phytochrome B (phyB) is clearly the most important photoreceptor 
in the vast majority of responses to shade, in some cases 
redundantiy with other members of its clade. In Arabidopsis, phyB 
mutants display a constitutive shade-avoiding phenotype that is 
characterized by long hypocotyls and petioles, reduced chlorophyll 
content, early flowering [48]. In the previous study, transgenic 
Arabidopsis over-expressing AtFPFl was also deemed to share the 
similar phenotype to phjB mutant, and the authors deduced that 
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the lack of phytochrome B leads to an enhanced responsiveness to 
GA [10]. Also they found that constitutively -FPFi-expression 
plants contained slightly higher amounts of GA4 and GA20 than 
wild-type plants did. Here the transcripts of key genes in 
gibberellin biosynthesis GA3oxi and GA20oxi were analyzed and 
proved to show litde diflFerence between wild and transgenic 
plants. The expression of PHTB was checked to have only 44% 
amount of that in wild- type. In addition, the transgenic Arabidopsis 
presented typical shade avoidance responses such as early 
flowering, longer hypocotyls and petioles, reduced chlorophyll 
content under high R: FR ratio light conditions. Though 
functional complementation assay should be developed further, 
suppressed expression of PHTB in transgenic Arabidopsis of over- 
expressing of GhFPFl may be the reason why shade avoidance 
syndrome (SAS) was induced presumably. 

In summary', identification of FPFl homologous genes in G. 
raimondii L. (D- genome) and G. arboreum L. (A- genome), also 
further comparing the sequences provided the groundwork 
necessary for understanding the distinctions and similarities 
between the sequences in the same genus of the different species. 
Strongly activated API and suppressed FLC expression in 
transgenic plants suggested some points for understanding the 
mechanism of GhFPFl to promote flowering. It seemed that 
transgenic Arabidopsis exhibited typical shade avoidance responses 
which might be caused by reduced expression of PHTB. 
Nevertheless, other proofs should be complemented in the future. 

Supporting Information 

Figure SI Pair-wise alignment of FPFl homologs 
DF10009151 (G. raimondii L.) and Garb22630 (G. 
arboreum L.). 
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(TIF) 

Figure S2 Paraffin section analysis of flower bud 
differentiation in G. hirsutum L. TM-1. A, B, C and D 

corespond to four developmental stages in shoot apices when there 
were two cotyledons, one, two, and three true leaves flattened. 
Images of representative paraffin sections were shown here. Yellow 
arrows pointed out tht; shoot apical meristem (SAM) of every 
stage. White and red arrows indicated vegetative bud primordium 
(VP) and floral bud primordium (FP) respectively. 
(TIF) 
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